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Abstract. We present recent results from spectroscopic data and modeling of the biconical ionized gas in the Seyfert-2 galaxy 
NGC 4388. A field of ~ 2.6 x 2.4 kpc centered on the nucleus has been observed by means of the modern technique of integral 
field spectroscopy. The analysis of more than two hundred spectra allowed to study the physical characteristics of the gas 
in the surroundings of the active nucleus. The South-West ionization cone, revealed by the [OIII] /15007/Hy8 excitation map, 
shows high emission line ratios not completely supported by simple photoionization. Composite models which account for the 
combined effects of photoionization and shock show that such high [OIII]/HyS line ratios are emitted by low density (n =30 
cm" 3 ) gas inside large (D>1 pc) shocked clouds (V s = 100 km s"" 1 ) reached by a relatively low flux from the active nucleus. The 
data of the VEELR in the North-East cone by Yoshida et al. (2002) have been modeled. The results confirm that photoionization 
is the prevailing mechanism, but nontheless weak shocks are under way between colliding clouds with small (< 1 pc) sizes and 
low densities no<100 cm"" 3 , moving outward at relatively low velocities (V s = 100 kms~'). 
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1. Introduction 

The prevailing picture of an AGN structure, the so-called 
Unified Model, provides for a central supermassive Black- 
Hole (BH), whose gravitational potential energy is the ultimate 
source of the AGN luminosity. Matter falling into the BH loses 
angular momentum through viscous and turbulent processes 
in an accretion disk, resulting into emission of photons in 
the ultraviolet, soft-X and hard-X ray wavelength domains 
( |Urry & Padovani 1995) . 

According to the Unified Model the central source is sur- 
rounded by a dense obscuring thick torus which is able to 
restrict the emergent radiation from the nucleus to a bipolar 
cone with an opening angle determined by the torus geometry 
(Pogge 1989 1. Imaging in optical emission lines shows in 
some cases high ionization regions with a conical and/or 
biconical morphology. These so-called ionization cones 
are direct evidence that radiation escapes anisotropically 
( |Mulchaey et al. 1994| i. 

In this frame we have carried out an observational campaign 
to study the circumnuclear and extranuclear environment 
of a spectroscopically selected sample of Seyfert galaxies 



( Rafa nelli et al. 19951 1 by means of the integral field spec- 
troscopy. This technique offers much more advantages than 
the traditional long-slit spectroscopy, because it allows us 
to observe extended regions of nearby galaxies providing 
simultaneously recorded spectral information (A) along two 
spatial directions (a and 6). 

Here we present results from the observations of the Seyfert-2 
galaxy NGC 4388, a nearby (z=0.0084) highly inclined spiral 
galaxy (i ~ 15°), believed located in the core of the Virgo 
cluster. Indeed observed distortions in the outermost optical 
isophotes are a clue of the possible interaction between this 
galaxy and the other members of the cluster. 
The IUE spectra of NGC 4388 revealed the presence 
of an UV source occulted from direct view and pro- 
ducing an ionizing continuum with a steep spectral 
slope (IFerland & Osterbrock 1986] |Kinney et al. 1991[ 
|Kinney et al. 1993) . 

Narrow-band Ha +[N II] and [O III] /15007 images pointed out 
the presence of emission line regions with highly ionized gas 
distributed into two opposite cones: a S-W cone extended to 
30" and a N-E plume extended up to 50" above the galactic 
disk. Spectroscopic analysis confirmed the photoionization by 
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Fig. 1. The position of the spectrograph over the galaxy. The x symbol indicates the assumed position of the central engine. 



the central power-law continuum as the principal excitation 
mechanism (jPogge 1988 Petitjea n & Durret 1993) . 

Kinematic data obtained with Fabry-Perot Interferometer 
dVeilleux et al. 1999b> revealed non-rotational blueshifted ve- 
locities in the extraplanar gas northeast of the nucleus, likely 
produced by a bipolar outflow. 

Using deep narrow-band imaging obtained with the Suprime- 
Cam mounted at the Subaru Telescope, Yoshida et al. (2002; 
hereafter YOS02) have recently found a very large emission 
line region (VEELR) extended up to ~ 35 kpc northeast of the 
nucleus. It consists of many gas clouds or filaments, part of 
which (within 12 kpc) are clearly ionized by the nuclear radi- 
ation. They claim that the most plausible origin of this ionized 
gas is tidal debris due a past interaction with a gas-rich dwarf 
galaxy. 

Moreover at radio wavelengths NGC 4388 shows a two sided, 
asymmetric structure with a flat spectrum central component, 
offset by Aa = 1 .2" and A6- 2.6" from the apparent optical 
nucleus, a diffuse blob to the North and an elongated feature to 
the South with a compact blob in it (Hummel & Saikia T9 9 1 1. 
Several radio maps have shown that this feature is clearly asso- 
ciated with the S-W ionization cone, suggesting the possibility 
of an interaction between the radio jet and clouds of the ion- 



ized gas (Falc ke et al. 19 98 1. Such interaction generates shock 
fronts, which cannot be neglected in the interpretation of the 
spectra. 

In this paper we investigate the physical conditions in the 
extended biconical ionized gas of NGC 4388 through the mod- 
eling of the optical emission line and continuum spectra ob- 
served in different regions. Both the ionizing radiation from the 
active nucleus and the shock effects are consistently accounted 
for. In Sect. 2 the observations and modeling of the S-W cone 
are presented. In Sect. 3 the N-E cone is considered and the 
VEELR is modeled on the basis of YOS02 data. The contin- 
uum SED is discussed in Sect. 4. Concluding remarks follow 
in Sect. 5. 

2. The South-West cone 

2.1. Observations and data reduction 

NGC 4388 was observed in March 1998 at the 6-m telescope 
of the Special Astrophysical Observatory (Russia) with the 
MultiPupil Fiber Spectrograph (MPFS), an integral field unit 
made by an array of 16 x 15 microlenses coupled with a close- 
packed bundle of optical fibers, which carry the incoming sig- 
nal to the spectrograph. Each of the 240 elements was covering 
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a region of l"xl" corresponding to a total field of view of 
about 2.6x2.4 kpc (assuming Ho = 75 km s _1 Mpc~'). Two 
600 sec exposures were taken using a TK-1024 CCD, which 
has a pixel size of 24 //m, and orienting the largest side of the 
array at P. A. ~34° (Fig.^. The grating was chosen to give a 
spectral range of ~ 4400-7 100 A, with a resolution of 10 A and 
a dispersion of 5.3 A/px. 

The reduction of the data was carried out by means of a 
software package, developed within our group, and running un- 
der IRAF 1 environment. 

After the usual steps of bias and dark subtraction, all the frames 
were cross-correlated, in order to compensate for the shifts be- 
tween them caused by mechanical flexions of the instrument. 
Then a flat field was chosen as reference image to obtain the 
geometrical scheme of the spectra arrangement: the position, 
width and layout along the dispersion direction of every spec- 
trum stored in it were defined and the resulting scheme was 
subsequently applied to all the images in order to extract for 
each of them the 240 1-D spectra. The dispersion solution for 
wavelength calibration was determined by fitting the line posi- 
tions of He-Ne-Ar lamp exposures with low order Chebyshev 
polynomial functions. The sky subtraction was performed by 
averaging the output of 8 fibers devoted to the night sky obser- 
vation and placed at a distance of 4.5' from the center of the 
field of view. Then the data were flux calibrated by observing 
the spectrophotometric standard star HZ 44: all its fiat-fielded, 
wavelength calibrated and sky subtracted spectra were summed 
together in order to collect the total incoming flux, to obtain 
the function necessary to convert counts in physical units and 
to account for the sensitivity response of the system at different 
wavelengths. Finally, cosmic rays were removed by combining 
corresponding spectra of the two exposures of NGC 4388 and 
a correction for Galactic extinction (Ay = 0.08) was applied. 
Fig-Elshows some examples of the reduced spectra chosen with 
different signal-to-noise ratio in the field of view. 
The positions, widths and fluxes of all the emission lines 
detectable 3<x over the continuum were measured using the 
IRAF task SPLOT. In case of close blending like Ha +[NII] 
ii6548,6583 or [S II] ^6717,6731, better results in perform- 
ing an interactive multigaussian fit of the total profile could be 
obtained by using the MIDAS package ALICE, whose input 
parameters (a first estimate of amplitude, center and sigma of 
the fitting gaussian) are independently fixable. 

2.2. Data analysis 

In order to get accurate estimate of line fluxes and ratios it is 
necessary to take into account that emission lines in galaxies 
are generally affected, and in some cases even dominated, by 
the underlying starlight contribution. It is possibile to remove 
the stellar features by subtracting an absorption-line "tem- 
plate" galaxy spectrum built by means of spectral synthesis 
techniques (see e.g. |Bonatto et al. (1989)) , or directly obtained 




5000 5500 6000 6500 7000 

Wavelength (A) 

Fig. 2. Examples of three 1-D spectra after the reduction. They 
are representative of different signal-to-noise ratio in the field 
of view. 




1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., 
the National Science Foundation 



under cooperative agreement with 
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Fig. 3. Example of the underlying starlight correction. An 
absorption-line "template" spectrum, conveniently rescaled, is 
subtracted from one of the 240 observed spectra in NGC 4388. 
The result shows that metallic lines are canceled out and the 
underlying H/3 absorption is removed. 



from the observation of early-type spirals. In our case the sec- 
ond approach was used, chosing a spectrum of a SO galaxy ex- 
tracted from our database which best matched the stellar fea- 
tures of most of NGC 4388 spectra, and following the sug- 
gestions given by Ho et al. (1993) As a whole, the resulting 
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correction was considered satisfactory (an example is given in 
Fig.gJ. The metallic lines (Mg I A5175, Fe 1 ,15269 and Na I 
/15892) were almost always canceled out after the subtraction. 
Only in some cases the stellar features of the template have 
been strengthened to better match the H/3 absorption, produc- 
ing an overcorrection of the metallic absorptions. This can hap- 
pen since stellar populations cannot be omogeneus everywhere 
in the field of view. The resulting variations of line intensities 
have been estimated weak (< 10% for H/3 ) in spectra with very 
bright emission lines, significant (~ 10 - 30%) in most of the 
spectra and very important in few spectra (up to 60%). 

The emission line fluxes, measured after the starlight 
contamination removal, were corrected for internal redden- 
ing by means of the theoretical Balmer decrement value 
Ha fH/3 = 2.86 for Case B recombination at electronic 
temperature T e =10 4 K and assuming the reddening law by 
IMiller & Mathews (1972)| 

Several images of the galaxy within the field of view were 
reconstructed at different wavelengths, as for instance the map 
of the visible continuum at 5500 A. The position of its maxi- 
mum intensity is located to the N-W of the observed peak of 
the emission line maps, which we have decided to assume as 
position of the photoionization source (hereafter nucleus). We 
do not confirm the large offset between the location of the radio 
source with flat spectrum and the visible continuum peak indi- 
cated in Hummel & Saikia (1991) and references therein: the 
radio source, the visible continuum and the emission lines have 
their maxima intensities close to each other (separation ~ 1"). 
We have discarded the idea to adopt the position of the radio 
source, which should indicate the real location of the AGN, as 
nucleus in our field of view since the spatial resolution of our 
data is only 1 "/px and an accurate astrometric calibration of 
the integral field images is very hard to obtain. 

The map of the [O III] /15007/H/? emission line ratios re- 
vealed part of the S-W ionization cone (Fig.0}. The cone has 
the axis at P.A. ~200° and an aperture angle of ~70° , less 
than the value measured by Pogge (1988), probably due to our 
smaller field of view. 

The gaseous regions within the ionization cone have been 
investigated. Their intrinsic FWHM, estimated from the gaus- 
sian analysis of the emission lines after having removed 
the instrumental component (~10 A), range in the interval 
~200-700 kms~', with a median value of ~480 kms~', ve- 
locities usually observed in the Narrow Line Regions of 
Seyfert galaxies. The diagnostic emission line ratios [Nil] 
^6583/Ha , [OI] A6300fHa , [S II] AA671 6 +673 1/Ha and 
[O III] ^5007/HjS JVeilleux & Osterbrock 19871 . show the typ- 
ical properties of gas ionized by a non-thermal source (Fig. [5}- 

Moreover the spatial distribution of the ionization degree 
within the cone, evaluated by means of the [O III] /15007/H/3 ra- 
tio as a function of the distance from the nucleus, indicates that 
the gas is highly ionized even at large distances from the AGN. 
In order to understand if the AGN could account for such ob- 
served ionization, the Ha luminosity of each spectrum was cal- 
culated from the reddening corrected Ha fluxes. Then the num- 
ber of ionizing photons needed to produce such luminosities, 
Q ion = 7.3 x 10 11 L(Ha) photons sec 1 dKennicutt T9 98 1, was 
evaluated and later compared with Q* uc = Q nuc Q/4w, where 



Q nllc is the number of photons coming from the nucleus, di- 
luted by the covering factor Q/47T, and Q. is the solid angle of a 
region seen from the nucleus. 

By assuming that the nuclear source is an isotropic emitter 
of radiation occulted along our line of sight, we attempted to 
estimate Q nuc by taking the average of the calculated number 
of photons ionizing the regions surrounding the nucleus, after 
having removed the effect of the covering factor (Qj on 4-n/Q). 
A value of Q nuc ~ 4.87 x 10 52 photons sec -1 was obtained, in 
good agreement with the 4.98xl0 52 photons sec -1 computed 
by Colina (1992) using the IUE spectrum of the galaxy and the 
flux at 912 A. 

The resulting Qion/Q* U c rat i° concentrates around a value of 
2.5. This indicates that the active nucleus is the dominant ioniz- 
ing source for the regions within our field of view, in agreement 
with what observed in the VO diagrams. In spite of that, the 
ionization parameter U = Q nU c/47rcrgNH, calculated within the 
cone as a function of the distance ro from the nucleus, ranges 
between 10~ 5 and 10~ 2 8 , a result clearly in contrast with the 
already pointed out high ionization level of almost all the re- 
gions. A total hydrogen density Nh = 3 x 10 2 cirT 3 is assumed, 
on the basis of the the electron density values estimated over 
the whole field, using the [S II] AA67 17/6731 ratio and assum- 
ing T=10 4 K (except where the sulphur lines were not visible or 
extremely blended). The problem of the missing data has been 
solved taking the average value of independent linear interpo- 
lations of the [S II] ratios along the two orthogonal directions 
of the spectral array. The resulting N e distribution presented 
values <10 3 cirT 3 , with a median of ~300 cirT 3 . Of course due 
to the uncertainties of the [S II] emission line strengths, the es- 



Fig. 4. The ionization cone obtained by the reconstructed map 
of the [OIII]/Hy6 emission line ratios. The x symbol indicates 
the position of the assumed nucleus (as in Fig. 1), OC is the 
maximum intensity of the optical continuum at 5500 A, and 
RS is the approximate location of the radio source with flat 
spectrum. 
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log([0 I>6300/Ha) log([N I]]A6563/Ho) log([S 1I]A6716 + 31/Hu) 

Fig. 5. Diagnostic diagrams of the emission line regions inside 
the ionization cone. The typical uncertainties of the ratios are 
indicated by the error bars in the bottom-right corner of each 
plot. 

timation of the electron density is not very accurate and just a 
simple decrease by a factor of 10 can account for the ionization 
parameter in the regions far from the nucleus. 

2.3. The models 

The contribution of the shocks to the line intensities is 
tested by comparing calculated with observed spectra. The 
code SUMA (Contini & Viegas 2001, and references thereini, 
which accounts consistently for the photoionization flux from 
an external source and for the shocks is adopted. 

The general model assumes matter bound gaseous and 
dusty clouds, which move outwards from the source and emit 
the line and continuum spectra. A shock front forms on the 
outer edge of the clouds, while the inner edge is reached by the 
radiation flux from the source. 

The input parameters are : the radiation flux intensity, 
Fh (in photons cirT 2 s _1 eV~' at 1 Ryd), the spectral index, 
a (= 1.5 for all models), the shock velocity, V s , the preshock 
density, no, the preshock magnetic field, Bo (= 10~ 4 gauss), the 
dust-to-gas ratio, d/g, and the geometrical thickness, D, of the 
emitting clouds. Cosmic relative abundances (Allen 1973} are 
adopted. 

Composite models lead to a complex distribution of the 
temperature and the density throughout the clouds. In radiation 
dominated (RD) models the effect of the flux from the active 
center (AC) prevails on the shock, while in shock dominated 
(SD) models the flux is absent (Fh= 0). When shocks and pho- 
toionizing fluxes act on the opposite sides of the cloud, two 
regions of relative high temperature appear near the edges (see 
Sect. 2.4). Higher temperatures correspond to higher shock ve- 
locities (T oc V s 2 ), while the radiation flux cannot heat the gas 
to > 2 - 3 10 4 K. This temperature leads to high +2 /0 oxygen 
fractional abundance, so the [OIII] line strength depends on Fh- 

The cooling rate (oc n 2 , where n is the density of the gas) 
increases with compression downstream of the shock leading to 
cool gas in the internal region of the cloud. The larger the cool 
region the strongest the low ionization level lines. The relative 
sizes of the hot and cool gas regions depend, therefore, on the 
shock velocity, on the intensity of the flux from the AC, on 



the preshock density, and on the geometrical thickness of the 
clouds. 

A large grid of models is calculated in order to select 
the most suitable parameters to the observational evidences 
for each region in the galaxy. Recall that a first choice is 
roughly obtained considering that the [OIII] A5QQ7 line flux 
depends on the intensity of the nuclear radiation flux, the [S II] 
AA61 17/6731 line ratio depends on the density, as seen before, 
and the intensity ratios of high to low ionization lines depend 
on the geometrical thickness of the clouds and on the shock ve- 
locity (see Contini & Viegas 2001). The input parameters are 
then changed consistently in order to obtain the best fit of all 
the observed data. 

2.4. [01 1 \]/Hp and [O \]/H/3 

The observations cover the SW cone providing line ra- 
tios in many different regions. We focus on [OIII] 
,U4959+5007/H/3 and [OI] ,U6300+6363/H/3 line ratios, be- 
cause [O III] lines are generally strong in Seyfert 2 galaxies and 
are emitted by gas ionized by the radiation flux from the active 
nucleus, while [O I] lines are emitted from gas at lower temper- 
atures and are generally strong in the presence of shocks. 
Moreover, Veilleu xet al. (1999a)| claim that dynamical pro- 
cesses such as entrainment by AGN-powered radio jets gen- 
erally have a stronger effect on the kinematics of the highly 
ionized [O III] emitting gas than on those of the low-ionization 
Ha emitting material. 

Therefore, by modeling the observed spectra it is possible 
to determine the distribution of the radiation flux intensity and 
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Fig. 6. Contour map of the visible continuum at 5500 A. The 
isointensities correspond to the values: 1.2, 1.5, 1.8, 2.2, 2.5, 
2.8, 3.0, 3.15, 3.25 (x l.e-15 erg cirT 2 s" 1 ). Filled circles corre- 
spond to the values of the observed [O III]ii4959+5007/H/3 ra- 
tios. The small "squared window" located at (6,9) indicates the 
assumed position of the ionization source. The map is orien- 
tated like in Fig. 4. 
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Table 1. The models in Fig. 7. 
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Fig. 7. The comparison of model calculation with ob- 
servations (filled triangles) for [OIII],Li4959+5007/H J S vs 
[OI]^6300+6363/Hy8 . See text for a detailed description. 



of other parameters, e.g. the shock velocity of the gas, the den- 
sity, and the geometrical thickness of the clouds in the observed 
region. The reliability of the models is checked comparing the 
observed spectral energy distribution (SED) of the continuum 
with model results (see Sect. 4). 

The observed [O III]/H/3 line ratios overlaid to the contin- 
uum of the galaxy at 5500 A are shown in Fig. 6 in a syn- 
thetic representation. The [O III]/H/3 distribution is not smooth. 
Particularly, high [O III]/H/3 ratios appear in some regions far 
from the active center. This suggest that the clouds in those re- 
gions have abnormal intrinsic conditions. 

In Fig. 7 [OIII]/Hy6 versus [OI]/H/3 from the observations 
(black triangles) is compared with model calculations. All the 
models correspond to V s =100 kms -1 . Higher velocities are 
less fitting, and velocities lower by a factor < 2 lead to similar 
results. The thick solid line refers to shock dominated models 
with line ratios increasing with decreasing D and no. The bulk 
of the data, however, is scattered throughout a region of the 
diagram better fitted by radiation dominated models. The thin 
solid line determines roughly the right edge of this region. The 
numbers in Fig. 7 refer to Tabled where the models are de- 
scribed. The dashed lines connect models calculated with dif- 
ferent D (increasing from left to right) but all corresponding to 
the same flux intensity Fh, while the dotted lines connect mod- 
els calculated with different F H (increasing from left to right) 
and corresponding to the same D. Fig. 7 shows that the trend 
of both types of models is different from the observed one, i.e. 
[O III]/Hj0 increasing with [O l]/H/3 . 

The increasing trend of the line ratios can be explained ei- 
ther by a slightly larger D and stronger Fh (models 1 in Table 1) 
or by very large clouds (D > 1 pc, models 2 in Table 1), corre- 
sponding to a lower density, and a flux intensity not exceeding 
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Fig. 8. The grid of the Fh values (filled pentagons) overlaid to 
the contour map of the ionization cone. The isointensities cor- 
respond to the [Om]/U4959+5007/H/3 ratios: 9.5, 10.5, 12.5, 
14, 15.5, 17, 19. The map is oriented as in Fig. 4. The square is 
like in Fig. 6. 



the characteristic values calculated for most of the observed re- 
gions inside the cone. The other observed lines at each position 
are too few and too uncertain to constrain the models, there- 
fore, we suggest that the clouds reached by a strong nuclear 
flux (model 1) are located in the central region of the cone, 
while clouds showing the same high [OIII]/H/? ratios, but char- 
acterized by a low density and a large geometrical thickness 
(model 2), are located in some external regions of the cone. 
Extended clouds of rarefied gas appear in the terminal region 
of jets and winds. 

The distribution of Fh in the different regions of the cone, 
as it results from modeling is shown in Fig. 8. The range of the 
corresponding ionization parameters U is 2 10~ 5 - 5 10~ 4 . 

To better understand the results obtained above, the distri- 
butions of the electron temperature (dotted line) and density 
(solid line) throughout a cloud are given for different models in 
Figs. 9 and 10. Each figure is divided by a solid vertical line in 
two halves. The horizontal axis scale is logarithmic and sym- 
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metric to have a comparable view of the two sides of the cloud. 
The shock front is on the left and the radiation flux from the 
source reaches the right edge. 
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Fig. 9. The case of a small D and relatively low Fh- The distri- 
bution of the temperatures and densities (upper panels) and of 
the most significant ion fractional abundances (lower pannels) 
throughout the clouds for different models (see text). 
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Fig. 10. The same as for Fig.|9]in the case of a large D and low 
n . 

The geometrical thickness D is a crucial parameter to low- 
level and, particularly, to neutral lines, because they are emitted 
from rather cold gas. Therefore, in the bottom pannels the dis- 



tribution of the fractional abundance of O +0 /O (dotted line), 
+1 /0 (short-dash line), +2 /0 (long-dash line), and H +1 /H 
(solid line) have been plotted. Fig. 9 corresponds to D= 3 10 16 
cm (0.01 pc) and Fig. 10 to D= 4 10 18 cm (1.3 pc). It can be 
noticed that two zones of gas emit the [O III] lines, one from 
the side of the shock, where collisional ionization prevails, and 
the other from gas ionized by radiation. The ionization rate due 
to the photoionization flux from the nuclear source, dominates 
in a large region of the cloud, particularly in models calculated 
with V s =100 kms -1 , where the gas density and optical thick- 
ness are rather low. A large region of cool gas appears inside 
clouds corresponding to a large geometrical thickness and low 
density, leading to relatively strong [OI] lines. 

2.5. [SUJ/Hp 

Fig. 11 shows [S II] Ai6717+6731/H/3 versus [S II] 
AA61 17/6731. The low ionization lines are strongly af- 
fected by shocks. For models accounting for the shock a 
large range of electron densities follows from compression 
in the downstream region. The lines, calculated integrating 
through the geometrical thickness of the clouds, account for 
the different conditions of the gas (see section 2.4). 

Some significant models, which are the same as those 
adopted to explain the [OHI]/H/3 ratios in Fig. 7, are indicated 
by open circles. They show a good agreement to the observed 
[SIT] /M6717/6731 line ratio between 1 and 1.3 in Fig. 11. 
The low density model gives [S II] /U67 17/6731 < 1.4. Notice, 
however, that the [S II] -U67 17+673 l/H/3 ratios are higher by 
a factor < 2 than the data for most models, indicating that the 
S/H relative abundance is lower than cosmic. This is a common 
characteristic to Seyfert 2 galaxies where sulphur is depleted 
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from the gaseous phase. The models are calculated adopting 
d/g (by number) = 10~ 14 - 1CT 13 , indicating that dust is indeed 
present in the NGC 4388 cone region. 

Reminding that compression increases with V s , models cal- 
culated with higher densities and velocities (V s =300 kms~') 
fit the data corresponding to low [S II] /M6717/6731 (<1), al- 
though their observed errors are large. Finally, models with 
V s =500 kms -1 hardly fit any data, suggesting that high ve- 
locity clouds are relatively few (see Sect. 4). 

The results show that the bulk of the shock velocities is 
lower than the velocities corresponding to the FWHM of the 
emission lines (200-700 km s ), indicating that strong shocks 
do not form in the observed region of NGC 4388, but rather 
low shock velocities V s = 1 00 km s~ 1 and low densities no< 1 00 
cm -3 dominate. 

3. The North-East cone 

We now consider the N-E cone and the very extended emission- 
line region (VEELR) discovered by YOS02 in order to have 
a larger view of the ionization conditions in NGC 4388. The 
nuclear region of the galaxy has been successfully investigated 
by Petitjean & Durret (1993), but only a few emission lines 
were observed at 21", 25", and 38" from the nucleus. 

YOS02 observed the [OIII] ,15007 and Ha +[NII] 
,1,16548,6583 fluxes from different clouds in the VEELR, 
which extends northeastward of the galaxy, very far outside the 
galaxy disk. The ionization cone is traced by the N-E plume 
and a high ionization cloud located at 10"- 15" northwest of 
the nucleus. They claim that there is a sudden change in the 
excitation state of the gas around the line at RA. 65°, which 
coincides with the extrapolation of the northern edge line of 
the S-W cone. Notice that these authors assume as distance of 
NGC 4388, the distance of the Virgo cluster (16.7 Mpc), which 
is different from that one used by us for the S-W cone (33.6 
Mpc). Therefore, all the distances should be scaled by a factor 
of 2. 

In Figs. 12, 13, and 14 we compare the results of models 
calculated with the code SUMA with the YOS02 data. The 
errorbars refer to the typical uncertainty of 20 % estimated 
by YOS02. In all the diagrams the fluxes of [OIII] versus 
Ha +[NII] (in erg cirT 2 s -1 ) are presented. Recall that models 
are calculated at the emitting nebula, while data are observed at 
Earth. Therefore, we have multiplied the data by d 2 /r 2 , where r 
is the distance of each cloud from the center of the galaxy and 
d is the distance of NGC 4388 from Earth. 

All models are calculated adopting a geometrical thickness, 
D = 3 pc, and Bo = 10~ 5 gauss, which is similar to the magnetic 
field in the ISM (B = 10 6 -10 5 gauss). 

Shock dominated models, calculated considering only heat- 
ing and ionization of the clouds by the shock (Fh=0), are pre- 
sented in Fig. 12. The models are labeled by the shock veloc- 
ity in kms~'. Large squares correspond to no=10 cirT 3 , small 
squares to no=2 cirT 3 . Relatively low preshock densities lead to 
the best agreement with the data, however, the SD models give 
a poor fit to the observed trend. This confirms YOS02 sugges- 
tion that gas excitation depends also on the flux from the active 
nucleus even at such large distances. We have used, therefore, 
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Fig. 12. The comparison of model results (corresponding to 
F H =0.) with the observations for [OIII] ,15007 (erg cirT 2 s -1 ) 
versus [NIIJ+Ha' (erg cirT 2 s -1 ) in the VEELR. Symbols are 
given in the text. 
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Fig. 13. The same as for Fig.^]for composite models calcu- 
lated with no = 10 cm -3 . 



RD models to fit the data. Recall that RD models are also com- 
posed because they account for the shock, even if the effect of 
the photoionizing flux prevails. 

We present in Figs. 13 and 14 radiation dominated models 
calculated with different velocities and fluxes, adopting no = 10 
cm -3 and 2 cm 3 , respectively. Higher densities do not fit. In 
Fig. 13 black triangles indicate models calculated with V s =50 
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Table 2. The best fitting models. 
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Fig. 15. The comparison of selected model results with 
the observations for [OIII] /15007 (erg cm -2 s -1 ) versus 
[NIIJ+Ha (erg cirT 2 s -1 ) in the VEELR (see text). 
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Notice that, roughly, all the data can be reached by the models 
represented by the grid. 

In Fig. 14 the models are calculated with V s =50 
kms -1 (black triangles), V s =70 km s -1 (small white triangles), 
and V s =100 km s -1 (large white triangles). The adopted fluxes 
range between log Fh= 7. and 9. It can be noticed that the cal- 
culated trends do not correspond to the observed one, however, 
some data are well reproduced by the models. In particular, the 
maximum and minimum observed [N II]+Ha are not explained 
by these models. 

We have chosen the best fitting models from Figs. 13 and 
14 and we present in Fig. 15 consistent modeling of the entire 
dataset. The selected models represented by open circles are 
described in Tabled 

The fluxes of [OIII] and Ho- +[NII] are given in Fig. 16 
as function of distance from the center. The data refer to the 
fluxes observed at the clouds. Fig. 16 shows that the line fluxes 
decrease with distance, however, different conditions coexist in 
the VEELR leading to the large scattering. 

The emission line flux from the cloud, identified by YOS02 
(see their Fig. 3) as C20, located at the extreme northern edge 
of the cone, is explained by a low flux (log Fh= 6.5) , no= 2 
cm" 3 , and a low shock velocity (35 kms -1 ). The largest Ha + 
[Nil] flux is emitted by the cloud identified as C26, located 
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kms , white triangles in the left side of the diagram corre- 
sponding to V s =30 kms^'and at the right to V s =70-90 km s -1 . 
The models are labeled by the values of log(FH). The maxi- 
mum [NII]+Ha value is explained by V s =150 kms -1 . Dotted 
lines join the results of models calculated with V s = 90 kms -1 . 
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Fig. 16. The fluxes of [OIII] (filled triangles) and 
[NII]+Ha (asterisks) observed at the clouds are plotted 
as function of the cloud distance from the center. 
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in the nebula at the southern edge of the eastern cloud. It cor- 
responds to relatively high Fh (log Fh=9.5) and high shock 
velocity (< 150 kms -1 ). Also cloud C25 is slightly dislocated 
from the correlating curve in Fig. 15. Therefore, it shows condi- 
tions different from those of the other clouds. Cloud CI shows 
the lowest distance from center and is modeled by log Fh= 9, 
V s = 70 kms" 1 , and no = 10 cm' 3 (Fig. II 1>. However, CI also 
corresponds to a shock dominated model (Eig.ll2> with V s =80 
km s _1 and no=10 cirT 3 . We have not enough data to constrain 
the model and we conclude that many different conditions co- 
exist in the clouds corresponding to the observations; in fact, 
the flux from the active nucleus, screened by the intervening 
matter, cannot reach all of the filaments within the observed 
clouds. The geometrical thickness adopted in the models is D= 
3 pc, while the observed clouds show dimensions of the order 
of 100 pc. Shock dominated models are radiation bound, so 
the results should not change even adopting lower geometrical 
thickness. 

To fit the data in Figs. 12, 13 and 14, we multiply the mod- 
els by factors of 10~ 5 , 10~ 8 and 10~ 65 , respectively. They ac- 
count for the filling factor. Cloud C20 has a surface area of S 
= 60x 190 pc 2 and the distance from center is r = 14.9 kpc 
JYoshida et al. 20021 Table 2). From S/(4 n r 2 ) x ff = 10~ 6 5 
the filling factor results ff = 0.08. For cloud C26, which cor- 
responds to gas with a higher density, we obtain ff = 1.3 10~ 4 . 
Regarding cloud C 1 , ff = 1 .4 1 0~ 4 results for the radiation domi- 
nated filaments and ff= 0.14 for the shock dominated filaments. 

YOS02 claim that the inner region (r< 12 kpc) of the 
VEELR may be excited by nuclear ionizing radiation, while 
the excitation mechanism of the outer region (r> 12 kpc) is 
unclear. Our analysis is not complete because the [O III] line 
fluxes were not observed for r> 15.6 kpc. It is found, however, 
that clouds C19-C22 correspond to lower fluxes and lower ve- 
locities, because the nuclear flux from the AC is diluted by dis- 
tance. The ionization parameter U that results from modeling 
the VEELR ranges between 10~ 5 and 1.5 10~ 3 . The shock ve- 
locities (50-150 km s~ 1 ) are in agreement with the radial veloci- 
ties in the N-E plume, i.e. 100-200 km s"' JYoshida et al. 2002> . 



4. The SED of the continuum 

The consistency of the models is checked in Fig. 17 where 
the SED of the bremsstrahlung radiation and re-radiation by 
dust from different clouds within NGC 4388 is compared with 
the data. Therefore, in Fig. 17 two lines correspond to each 
model, one indicating the bremsstrahlung from the gas and the 
other emission from dust in the infrared. The observational 
data (from the NED 2 ) are taken from Boroson et al. (1983), 
Scoville et al. (1983), de Vaucouleurs & Longo (1988), Soifer 
et al. (1989), Moshir et al. (1990), Becker et al. (1991), de 
Vaucouleurs et al. (1991), Gregory & Condon (1991), Kinney 
et al. (1993), Fabbiano et al. (1992), Zwicky et al. (1961), 
Douglas et al. (1996), Dressel & Condon (1978), Spinoglio et 
al. (1995) and Mould et al. (1980). The two data in the soft X- 
ray range come from Einstein (Fabbiano et al. 1992), integrated 

2 NASA/IPAC Extragalactic Database 



t— , ,— | r— i i— | i—i i— | i—i i— | i—i r 




log(i/) 

Fig. 17. The SED of the continuum. The data are represented by 
filled squares. V s = 100 km s -1 , n = 100 cm" 3 , logF H = 8.77, 
D=0.01 pc, d/g= 1.5 10~ 14 (short-dash lines), V. s =100 km s _1 , 
n =30 cm 4 , D=1.3 pc, log F H = 8.95 d/g=10~ 14 (dotted lines), 
and V s = 500 km s _1 , n = 300 cnT 3 , D=3.3 pc, d/g= 1.5 
10 -13 (long-dash lines). The summed continuum SED (thick 
solid lines); radio synchrotron continuum (thin solid line). 



on 0.2-4.0 keV, and from ROSAT (Halderson et al. 2001), inte- 
grated on 0.1-2.4keV. 

Indeed, the data from the NED were taken with differ- 
ent precisions and in different epochs. Nevertheless, they are 
very significant because they correspond to the main contribu- 
tions to the continuum in the different frequency ranges by 1) 
bremsstrahlung from the photoionized gas (v ~ 5 10 13 - 5 10 15 
Hz), 2) bremsstrahlung from the hot gas in the downstream re- 
gion of shocked clouds (v > 10 16 Hz), 3) emission from the old 
star population (v ~ 5 10 I3 -10 15 Hz), 4) dust reradiation in the 
infrared range (v ~ 10 12 -10 14 Hz), and 5) radio emission (v < 
10 11 Hz). 

Besides the line and continuum spectra, the SUMA code 
calculates reradiation from dust in a consistent way. Dust 
is heated by the nuclear radiation flux; moreover, gas and 
dust mutually heat each other by collisions in shock turbu- 
lent regimes (Viegas & Contini 1994). Therefore, the frequency 
corresponding to the dust reradiation peak depends on the 
shock velocity. At high V s the grains will reach high temper- 
atures in the post shock region and dust reradiation will peak 
in the mid-IR, while at low velocities the peak appears in the 
far-IR. The peak intensity depends on the dust-to-gas ratio. 

Sputtering of the grains, which depends on the shock veloc- 
ity, is calculated throughout the clouds. Small grains are rapidly 
sputtered, therefore we adopt silicate grains with an initial ra- 
dius of 0.2 fan. 
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Emission in the radio range corresponds to bremsstrahlung 
from cool gas as well as to synchrotron radiation due to Fermi 
mechanism at the shock front. 

The continuum SED corresponding to single clouds are 
shown in Fig. 17 in order to explain the role of the different 
models. Then, the contributions from all the clouds are summed 
up to fit the data consistently in all the frequency domains. 

The models obtained for the S-W cone are used to model 
the SED of NGC 4388, because the emission in the VEELR is 
very low and can hardly affect the continuum observed from 
the entire galaxy. The models presented in Fig. 17 are se- 
lected from those which better fit the bulk of the data shown 
in Fig. 7. We will consider three significant models, regarding 
the shock velocity, the preshock density, the geometrical thick- 
ness of the cloud, and the dust-to-gas ratio. The short-dashed 
lines refer to a model calculated with V s =100 kms" 1 , no=100 
kms 1 , D=0.17 pc and d/g=1.5 1(T 14 . Notice that d/g = 1(T 14 
by number corresponds to 4.1 10~ 4 by mass (Draine & Lee 
1984) which is the dust-to-gas ratio in the Galaxy adopting a 
silicate density of 3 g cirT 3 . Dotted lines represent the SED of 
the model calculated with D=10 pc, no=30 cm 3 , log Fh= 9.8, 
and d/g=l(T 14 . Considering that the FWHM of the line profiles 
show velocities up to > 500 km s _1 , we have added in Fig. 17 
also the continuum of a SD model corresponding to V s = 500 
kms -1 , no=300 cm 4 , and d/g =1.5 10~ 13 , in order to fit the 
data at higher frequencies. This relatively high shock velocity 
leads to a temperature in the immediate postshock region of ~ 
3.7 10 6 K. Consequently bremsstrahlung radiation peaks in the 
soft X-ray domain. 

Indeed, the analysis of the ROSAT HRI observations by 
Matt et al (1994) has shown that the 0.1-2.4 keV X-ray emis- 
sion is dominated by an extended component (~ 4.5 pc). It 
is interpreted as thermal emission from a collisionally heated 
plasma and a collection of discrete sources. The corresponding 
temperature of ~ 0.4 keV (~ 4.6 x 10 6 K) is in agreement with 
the temperature of ~ 4 10 6 K calculated by V s =500 kms" 1 . No 
correlation of the soft X-ray emission with the optical ioniza- 
tion cones has been found (Matt et al.). This confirms our pre- 
vious results (Sects. 2 and 3) showing that rather low (~ 100 
km s _1 ) shock velocities dominate in the cones. 

Fig. 17 shows that the SED of single models calculated 
with V s =100 kms -1 nicely fit the data corresponding to the 
lowest fluxes in the optical-UV range. The data between 10 14 
and 10 15 Hz are, however, contaminated by radiation from old 
population stars and are explained by a black body model cor- 
responding to 5 10 3 K (dash-dotted line). This is common to 
the continua of Seyfert 2 galaxies ( jContini & Vie gas 2000 1. 

Recall that dust and gas mutually heat each other across 
the shock front, therefore, dust re-radiation from high ve- 
locity clouds (V s = 500 kms -1 ) peaks in the mid-IR (see 
IContini et al. 20021 . Our results show that these clouds are 
rather dusty. 

The data in the radio range do not follow exactly the slope 
corresponding to bremsstrahlung from a cool gas. Yet, they do 
not fit a power-law distribution (thin solid line), which could 
be explained by synchrotron radiation created by Fermi mech- 
anism at the shock front. Fig. 17 shows that both mechanisms 
contribute to the radio emission flux. 



Finally, the thick solid lines in Fig. 17 represent 
the weighted sum of the single-cloud models for gas 
bremsstrahlung and dust reradiation, separately. The line is 
interrupted at frequencies lower than 10 10 Hz, because self- 
absorption, which increases with wavelength, strongly reduces 
bremsstrahlung emission in the radio domain. 

The model best fit to the data is obtained adopting differ- 
ent weights for different models. The weights account for the 
relative number of clouds corresponding to the model. 

The ratio of the weights used to fit the data by the three 
models is 1 : 0.2 : 0.03, the maximum weight corresponding 
to V s =100 km s~'and D=0.17 pc and the minimum to V s =500 
kms -1 . This confirms that velocities of ~ 100 kms^'and rela- 
tively small D (< 1 pc) dominate. 

5. Concluding remarks 

The central regions of the Seyfert-2 galaxy NGC 4388 have 
been observed by means of the integral field spectroscopy tech- 
nique. The [O III] /15007/HjS map revealed part of the S-W ion- 
ization cone, with P.A. ~200° and aperture angle of ~70°. 
The analysis of the spectra within the cone shows the non- 
thermal active nature of the emission line regions, whose ion- 
ization degree remain high and almost constant even at large 
distances from the nucleus. With simple energy balance cal- 
culations we estimate the number of the ionizing photons pro- 
duced by the engine and we demonstrate that even if the pho- 
toionization is the main physical process responsible for the 
excitation of the gas inside the cone, it cannot be considered 
the only one under way, since it cannot explain the observed 
high ionization far from the source. 

We verify our hypothesis by applying to our spectroscopic 
data the SUMA code, which accounts for the combined ef- 
fect of the photoionization from an external source and of the 
shocks. The results of the models, whose reliability is cross- 
checked by comparing them with the observed SED of the 
continuum, show that the ionization cone is characterized by 
a mixture of low density colliding clouds with large (D > 1 pc) 
and small (D<0.1 pc) geometrical thickness. Collisions, in fact, 
cause turbulence and R-T and K-H instabilities which lead to 
fragmentation of matter. In particular we are able to explain the 
high [OIII]/Hy6 line ratios observed in the outer edges of the 
S-W cone by the coupled effect of a low density gas in large 
clouds (D > 1 pc). 

The models show that the bulk of the emitting clouds within 
the S-W cone has shock velocity of ~ 100 kms 1 , definitively 
lower than those corresponding to the observed FWHM of the 
line profiles (200 - 700 kms 1 ). Adopting the hypothesis that 
the extraplanar material in NGC 4388 may represent the tidal 
debris from a recent encounter with a small gas-rich galaxy 
(Pog ge 1988|IVeilleux et al. 1999bt . we suggest that head-on- 
back shocks accompany the outwards motion of the clouds, 
which result when an incoming gas catches up with preexistent 
gaseous clouds moving in the same direction with a different 
velocity. 

To complete the analysis of the biconical structure of NGC 
4388 we have modeled the emission line fluxes observed by 
YOS02 in the VEELR of the N-E cone. We find that the phys- 
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ical conditions are in agreement with those found in the S-W 
cone, taking into consideration dilution with distance from the 
AGN. The results lead to log F H = 8.3 - 10.5 (U = 2 1(T 5 - 5 
1(T 4 ), n = 30 - 100 cm" 3 , and V s ~ 100 -300 kms" 1 in the 
S-W cone, and log F H = 6.5 - 9.5 (U = 10~ 5 - 1.5 10" 3 ), n = 
2 - 10 cm 3 , and V s = 30 - 150 kms 1 in the VEELR of the 
N-E cone. We prefer to compare the radiation fluxes instead of 
the ionization parameter, because U depends on the density of 
the gas in the edge of the cloud facing the radiation source. So, 
considering that the density decreases outward the galaxy, di- 
lution is less prominent. The values of U in the nuclear region, 
calculated by Petitjean & Durret (1993), are much higher, U = 
5 10~ 3 - 0.014, indicating that different conditions characterize 
the gas in this region. 

Concluding, densities, shock velocities, and the magnetic 
field are lower in the VEELR, suggesting that beyond the edges 
of the VEELR the emitting clouds are merging with the ISM. 
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